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Abstract 

I propose an orthogonalization procedure preserving the grading of the initial graded 
set of linearly independent vectors. In particular, this procedure is applicable for 
orthonormalization of any countable set of polynomials in several (finitely many) in- 
detcrminates. 

There are two well-known procedures for orthogonalization of any set of linearly indepen- 
dent vectors in a linear space. The first of them boils down, essentially, to calculation of 
the (-s)-th power of the Gram matrix of the initial set (we refer to it as Gram method), 
and another one is the Gram-Schmidt process. The result of the Gram-Schmidt process 
essentially depends of the order of the elements in the set and the operations with Gram 
matrix seem to be impossible to perform in infinite dimensional spaces. 

There are, however, some situations when one has to orthogonalize a finite or infinite 
set of linearly independent vectors in a linear space in such a way that the result would 
have some properties of the initial set. Observe that the set considered is not necessarily 
a basis. Any set consisting of linearly independent vectors will do. Observe also that I 
deal here with algebraic problems, no analytical problem (convergence, completeness, etc.) 
arises. 

We consider some examples: 

1. Let {e m = e zmx ^°°__ be a set of linearly independent vectors in L2QO, 2tt], p), 
where p is any positive weight. If p is not constant, the elements of this set are not 
orthogonal, but have an important property e m = e_ m , where the bar means complex 
conjugation. Is it possible to preserve this property under orthogonalization? 

2. Let x = (xi, . . . ,x n ), where Xj £ R, and m = (mi, . . . , m n ), where m,j G Z+. 
Then the set {x m : m £ Z™ }, where x m = x™ 1 • • • x™ n , is, due to Stone- Weierstr ass 
theorem, a basis in £2(0) for a "good" bounded domain C W l . This basis has 
no natural order, but has a natural grading: the degree of x m is equal to \m\ = 
m\ + • • • + m n . How to orthogonalize this basis and preserve the natural grading? 
Is this possible? 
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It is very strange that a very simple and natural answer to these and similar questions 
seems to be unknown. For this reason I present here the orthogonalization procedure, 
which combines some features of the Gram method and Gram-Schmidt process (in par- 
ticular, it can be applied in infinite dimensional case) and gives a solution to the above 
problems and similar ones. Note in this connection that, although the orthogonalization 
problem is more than hundred years old, various aspects of the propblem arise from time 
to time in connection with very interesting practical problems, see, e.g., Q. 

Let {e^, : a £ Ik, k £ Z + } be a set of linearly independent vectors in a Hilbert space. 
We assume that all index sets Ik are finite. Let us inductively define the sets of vectors 
{f k : a € Ik} by the formula 

f k a = E Qk4 + EE C 1 ) 

P€h 3=0 Peij 

where the unknown matrices Q k and P k i are determined from the orthonormality condi- 
tions for the system /: 



ft) = $al3, a,(3 e I k , 

(2) 

=0, a el k , pel,, j + k. 

Note, first of all, that the system / is linearly independent if and only if the matrices Q k 
are nondegenerate. We will use this fact later. 

From the last line in (^) and the definition (Q) one immediately deduces that 

P& = ~ E ( 3 ) 

where we define the matrix D kj to be D k J = (e k ,fi). Let us substitute this expression 



for P in terms of Q and D into the first line in (|2|). We obtain, after simplification, a 
matrix equation of the form 

k-l 

Qk] B kQk = E ^ where B k =T k ( 4 ) 

3=0 

and where f denotes the Hermitian conjugation, A fcj = D kj ^D k3 , T^a = (e^,e^) is the 
Gram matrix of the system {e k Y } a( zj k and E is the unit matrix. 

Note that the matrix B is the Gram matrix for the linearly independent system of 
vectors in the Hilbert space 

^ = 4-EE( e -4)* 

j=o Pel, 

Hence, B k is positive definite. So one can write the unique positive definite solution of 
the equation (|j) in the "Gram" form 

( k ~ l \ ~ 1 ' 2 

Q k = lr k -^2 Akj ■ (5) 
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This completes the ortho-normalization process. 
We consider now some simple examples. 

3. If for all k the sets Ik are one-element sets, then the process described is exactly the 
Gram-Schmidt one. 

4. In the above procedure Z + can be replaced with any its finite subset. If such subset 
is a one-element set, then our process is exactly the Gram method. 

These two examples show that the process suggested is simply a combination of two 
well-known processes. 

5. Let Iq = {0} and = {+,— } for nonzero k's. Let further e(* G ^([0, 2tt], p) and 
e° = 1, e± = e ±linx . (This is Example 1.) Then, clearly, /° = const G R. Let us 

show now that if /+ = f 3 _ for all j < k, then the same is true for k also. 
Indeed, due to relations @ and (||) we obtain 

f± = Qi + (4 - (4. /°) /°) - E ((4. A) A + (4. A) A) 

3=1 

+ Qi- (e k . - (etj°) f) - J2 fi) fi + (>_, fi) fi) ■ 

Now using inductive hypothesis and the fact that Q are Hermitean matrices, i.e., 

G R and Q+_ = Q- + , it is easy to see that the vectors are complex conjugate. 
So we have obtained a simple affirmative answer to the question in Example 1. 
It is interesting whether or not the above construction can be generalized to the vector 
systems in pseudo-euclidean spaces? Two problems arise in this case: 

1) What shall we do with "isotropic" vectors? 

2) For which matrix in the right hand side instead of the identity one, is equation (Q) 
solvable? 

Note in this connection that the Gram-Schmidt orthogonalization process is not ap- 
plicable in the pseudo-euclidean case because even if the initial vector system does not 
contain isotropic vectors, such vectors can appear under execution of the process and 
terminate it. 

We consider an example: if vector e\ is isotropic, there does not exist any number 
a such that vector e\ + ae2 is pseudo-orthogonal to e± (of course, unless &\ and e2 are 
initially orthogonal). 

Certain properties of bases and linear independent systems in pseudo-euclidean spaces 
are discussed in 0, §]. In particular, Bognar proposes a "forceful" method for general- 
ization of Gram-Schmidt process to pseudo-euclidean case for linearly ordered systems of 
vectors We use something like his method for graded systems. 

In what follows we suppose that the pseudo-euclidean scalar product is "nondegene- 
rate", i.e., for any linearly independent finite system of vectors consisting of more than 
one element its Gram matrix is nondegenerate. In other words, this means that the 
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dimension of any maximal isotropic subspace (all of whose vectors are isotropic) does not 
exceed 1. Then any finite system of vectors, excluding one-element ones, may be "pseudo- 
orthonormalized" in the sense that the pseudo-norm, or "length" , of each final vector will 
be equal to ±1. 

To this end, it suffices to determine the signature (p, q) of the Gram matrix T of the 
initial system of vectors and then solve the equation 

R^TR = E P ® (-E q ) (6) 

for unknown matrix R. If we write the Hermitean matrix T in the form T = UAU^ , where U 
is a unitary matrix and A = A p © (— A q ) is a diagonal matrix with p positive and q negative 

elements, then the solution of (|^) can be written in the form R = U © ^— A q • 

Therefore, it is clear that the only trouble which might appear when dealing with 
graded systems of vectors in pseudo-euclidean spaces is when after the k-th step the set 
{fa '■ a € Ik} consists of exactly one isotropic vector. There does not exists any general 
natural way to resolve this situation and details depend on the concrete case. A most 
simple idea is to include such a vector into the next, (k + l)-th, level set of vectors and 
then continue the process. This trick slightly violates the initial grading, but preserves 
the filtration. 

The problems discussed stem from several sources. One is numerical analysis and data 
processing. Here the role of discrete Fourier transformation is well-known, but what should 
we do if our data is given on a nonuniform grid? For the answer see Example 1. 

Another problem is orthonormalization of splines which constitute a not linearly or- 
dered, but a graded set of functions. 

D Leites and A Sergeev pointed out a totally different area in which the same question 
arises. These problems concern with new polynomials in several indeterminates connected 
with some Lie algebras and superalgebras and the space of these polynomilas is naturally 
endowed with a nondegenerate indefinite metric, see ||. For one indeterminate D Leites 
and A Sergeev can orthogonalize their polynomials; for several indeterminates these poly- 
nomials are not linearly ordered and they got stuck. By our method one can orthogonalize 
the polynomials in several indeterminates proposed in Q. 

For the reader who wishes to compare various orthogonalization methods I suggest very 
transparent and user friendly paper by Srivastava 
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